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1. Introduction 

Recent experiments at Brookhaven's Relativistic Heavy Ion Collider (RHIC), give strong evi- 
dence for the creation of locally thermalized hot and dense QCD matter in a "little bang" initiated 
in ultra relativistic nucleus nucleus collisions The energy density e, achieved a short time after 
the little bang (~ Ifm/c) has been estimated to be > 5.4 GeV/fm 3 [Jl|] — a value significantly larger 
than the ~ 1 GeV/fm 3 required for the transition from hadronic matter to the high-temperature 
plasma phase of QCD [g, g|. A confluence of experimental results [|l], |], |], |6|] now suggest that this 
new state of matter bears strong kinship to the hot and dense plasma of quarks and gluons (QGP) 
produced a few micro seconds after the "big bang" that gave birth to our universe some 12 — 14 
billion years ago. Now, the scientific challenge is to explore robust experimental constraints which 
can establish the detailed properties of this QGR 

In this contribution, we show experimental validation for; (i) the scaling predictions of perfect 
fluid hydrodynamics for the elliptic flow coefficient V2, (ii) the development of elliptic flow in the 
pre-hadronization phase, (iii) the scaling of D meson flow compatible with full thermalization of 
the charm quark and (iv) universal scaling of the flow of both mesons and baryons (over a broad 
transverse kinetic energy range) via quark number scaling. Subsequently, we use the scaled flow 
values to constrain a number of transport coefficients. 



2. Thermalization and harmonic flow correlations 

Local thermalization plays a central role in considerations involving perfect fluid hydrody- 
namics. Consequently, it is important to recap a few of the key observables which are consistent 
with the production of high energy density thermalized matter at RHIC [0]. In brief; 

• the multiplicity distribution for the same number of participants is independent of colliding 
system size as would be expected from a system which "forgets" how it is formed. A demon- 
stration of this independence hypothesis has been given for Au+Au and Cu+Cu collisions by 
the PHOBOS collaboration [§]. 

• A comparison of the measured anti-particle to particle yield ratios p/p, to the predictions of 
a statistical model 

p e -(E-a)lT e ' {ZA) 
show excellent agreement for a single temperature T and baryon chemical potential jj, 

• Surprisingly large harmonic flow correlations compatible with the development of sizable 
pressure gradients in a locally thermalized fluid are observed for a broad selection of hadrons 



comprised of light- and heavy quark combinations (tt, A^p, A,I2, 0,S,d,D, . . .) [1C, 11, 12, 



13, |14J. We emphasize here that rather strong interactions are required in this fluid to gener- 
ate significant flow for hadrons comprised of a heavy quark, ie. the heavy quark relaxation 
time is much longer than that for the light quarks. 

The magnitude of harmonic flow correlations is commonly characterized by the Fourier coefficients 



v„ as Q15J]: 

Vn = ( e K<P-*«)) = (cosnto-d)*)), (2.2) 
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where <p is the azimuthal angle of an emitted particle (measured in the laboratory coordinate sys- 
tem) and is an estimate of the azimuth of the reaction plane. A correction factor for v„ is usually 
required to take account of the dispersion of the estimated reaction plane [16, 17, 18]. The second 
harmonic coefficient V2 is termed elliptic flow. 

At RHIC energies, it is widely believed that elliptic flow results from hydrodynamic pressure 
gradients developed in an initial "almond-shaped" collision zone. That is, the initial transverse 
coordinate-space anisotropy of the collision zone is converted, via particle interactions, into an az- 
imuthal momentum-space anisotropy. Elliptic flow self-quenches due to expansion of the collision 
zone. Consequently, rapid local thermalization is required to achieve relatively large V2 signals. 
The current consensus is that the sizable values observed at RHIC are compatible with the hydro- 
dynamic expansion of a near perfect fluid. 



3. Universal scaling and perfect fluid hydrodynamics 

Perfect fluid hydrodynamics stipulates that the constituents of a fluid flow with the same ve- 



locity field. In many hydrodynamic models [|19|, [20[ [21|, |22J, |23|, [24|, [25j], elliptic flow results from 
pressure gradients due to the initial spatial asymmetry or eccentricity e = ((y 2 —x 2 ))/((y 2 +x 2 )), 
of the high energy density matter in the collision zone. The initial entropy density S(x,y), is often 
used to perform an average over the x and y coordinates of the matter in the plane perpendicular to 
the collision axis. For a system of transverse size R(l/R = ^\/(x 2 ) + \/(y 2 )), this flow develops 



over a time scale ~ R/ (c s ) for matter with sound speed c s [g5]. Thus, the initial energy density 
controls how much flow develops globally, while the detailed development of the flow correlation 
patterns are governed by e and c s . Centrality dependent differential measurements give clear ac- 
cess to e; access to c s is less direct but is reflected in the magnitude of the flow. Lattice QCD 
calculations indicate that c 2 ^ \/3 for T > 2T C , but drops steeply by more than a factor of six near 
the "softest point" where T T c ie. the critical temperature p6[]. Thereafter, it rises again in the 
hadron resonance gas phase to the value c 2 « 0.15 for SPS energies. 

A remarkable empirical fact recently reported, is the observation that universal scaling for 



the flow of both mesons and baryons is achieved when V2/n q e is plotted vs. KEj/fiq [Zq, 
Here, n q is the number of valence quarks (n q = 2, 3 for mesons and baryons respectively) and 
KEj = m t — m is the transverse kinetic energy 1 . Further demonstration of this scaling is illustrated 
via Figs, [j] and ||. Figure [j] shows differential flow measurements V2 (pr ) , for several particle species 
produced at mid-rapidity in central and semi-central Au+Au collisions at ^snn = 200 GeV; they 
span essentially the full range of measurements (several hundred data points) at RHIC. The V2 
values shown for charged and neutral pions (n ± , 7T°), charged kaons (K^), (anti-) protons ((p)p), 
(anti-) deuterons {{d)d and the <p meson represent recent measurements (published and preliminary 
results) obtained by the PHENIX collaboration [0, |27|, |j| ||]. The values for neutral kaons (K®), 
lambdas (A + A), cascades (S + S) and omegas (Q. + £2) show results from the STAR collaboration 



Hn|,|i3|, y,£oj]. 

The most striking feature of this plot is the large magnitude of V2{pr) observed for all particle 
species. Fig. || shows the scaled results obtained from the same data. It indicates that the relatively 



The pressure in an ideal gas can be seen to be a measure of its kinetic energy density. 
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complicated "fine structure" of V2 (ie. its detailed dependence on centrality, transverse momentum, 
particle type, etc) for particles produced at mid-rapidity can be scaled to a single function. 
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PHENIX( Phys.Rcv.Lett.91, Preliminary: QM0S, GRC 06 ) 

• - 7t + +7f : min.bias,0-10%,10-20%,20-30%,30-40%,40-50%,20-60% 
O - 71 • min.bias 

■ - K + +K": min.bias, 0-10%,10-20%,20-30%,30-40%,40-50%,20-60% 
*- p+p I min.bias, 0-10%,10-20%,20-30%,30-40%,40-50%,20-60% 

- (J : min.bias, 10-50% 
A - (j) : 20-60% 

STAR (Phys.Rev.Lett. 92, Phys.Rev.C72(2005), Preliminary QM05.SQM06) 

- 7T+7I" • min.bias 

□ - Kg : min.bias, 5-30%,30-70% 

^ - P+P ' min.bias 

- A+A : min.bias, 5-30%, 30-70% 

* - S+S : min.bias 
A - Q.+Q. : min.bias 



Figure 1: V2 vs. pj for several particle species produced at midrapidity in central and semi-central Au+Au 
collisions at Jltm = 200 GeV. The centrality and particle species selections are indicated. 




PHENIXf Phys.Rev.Lett.91, Preliminary: QM05, GRC 06 ) 

• - 7l + +7l": min.bias,0-10%,10-20%,20-30%,30-40%,40-50%,20-60% 

o- 71° : 

■ - K + +K": 
*-p+p: 

-d : 

A-<j) : 

STAR 



min.bias 

min.bias, 0- 10 % ,10-20 % ,20-30 % ,30-40 % ,40-50 % ,20-60 % 
min.bias, 0- 10 % ,10-20 % ,20-30 % ,30-40 % ,40-50 % ,20-60 % 
min.bias, 10-50% 
20-60% 



(Phys.Rev.Lett. 92, Phys.Rev.C72(2005), Preliminary QM05,SQM06) 



• min.bias 



□ - Kg : min.bias, 5-30%,30-70% 

Q - P+P * min.bias 

- A+A : min.bias, 5-30%, 30-70% 

* - S+S I min.bias 

A - Q.+Q, : min.bias 



0.5 1 
KE T /n q (GeV) 

Figure 2: V2/n q e vs. KEj /n q for several particle species produced at midrapidity in central and semi-central 
Au+Au collisions at Jsnn — 200 GeV. The data selections are indicated. 
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Figure 3: V2 vs. for charged hadrons obtained in Au+Au (a) and Cu+Cu (b) collisions for several 
centralities as indicated. Panel (c) shows Vi{centrality, pr) divided by k times the /^-integrated value 
V2(centrality), for Au+Au and Cu+Cu (k — 3.1 see text and Refs. [[281 Egh). 



On the one hand, this universal scaling validates the predictions of perfect fluid hydrodynamics 
[§§§,110]. That is; 



• V2(pr)/e should be independent of centrality; 

• V2(pr) should be independent of colliding system size for a given eccentricity; 

• for different particle species, V2(KEt) at mid-rapidity should scale with the transverse kinetic 
energy KEj = mj — m, where mj is the transverse mass of the particle; 

• v a {pt) x vKpt)- 



On the other hand, it corroborates the n q scaling expected from quark coalescence models Q35| , |3 



37| . [38p suggesting that quark-related degrees of freedom are present when elliptic flow develops. 

These separate aspects of the observed universal scaling are illustrated in Figs.[| and |]. Figures 
|5|a and ||b show the differential V2(pt) for charged hadrons obtained in Au+Au and Cu+Cu colli- 
sions respectively. The expected increase in V2(pr) as collisions become more peripheral and the 
Pt increases ifj, ||] is unmistakable. Fig. ||c show the resulting £-scaled values for both systems. 
The eccentricity is obtained from the pj -integrated V2 values i.e. e = k x v%. A Glauber model 
estimate of e in Au+Au collisions [ p 7 ^ ] gives k = 3.1 ±0.2 for the cuts employed in the analysis 
plph 



5 



What do elliptic flow measurements tell us about the matter created at RHIC? 



Roy A. Lacey 



0.3 



0.2 



0.1 



(a) 



~r 



IK 



T 



T 



• rr + +rr (PHENIX) Q p+p (PHENIX) 
■ K + +K" (PHENIX) O A+A (STAR) 
K° (STAR) □ ="+=* (STAR) 




| I I I I | I I I I | I I I I | 




1 2 3 

(GeV/c) 



1 2 3 
KE,. {GeV) 



I i i i i I i i i i I 



,05 



0.5 1 1.5 2 
KE/n, (GeV) 



Figure 4: v'2 vs. pj (left panel) and KEY (middle panel). The scaled results in the right panel is obtained 
by n q scaling the data shown in the middle panel. Results are shown for several particle species produced in 
minimum bias Au+Au collisions at Js~m = 200 GeV [E9fl. 



These scaled values shown in Fig. [|c validate the expected independence on colliding system 
size []7j, ^5]] and suggests that rapid local thermalization [40, |d|] is achieved. They also allow an 
estimate of c s because the magnitude of V2/£ depends on the sound speed c s [25]. A recent estimate 
[28, 29] gives c s ~ 0.35 ±0.05; a value which suggests an effective equation of state (EOS) which 
is softer than that for the high temperature QGP [26]. It however, does not reflect a strong first order 
phase transition in which c s = during an extended hadronization period. Such an EOS could very 
well be the reason why V2{pr) is observed to saturate in Au+Au collisions for the collision energy 
range = 60 - 200 GeV (4^1 . 

The other aspects 2 of the observed universal scaling are illustrated in Fig.[|. The left panel 
show measurements of the pj dependence of V2 for several particle species. The middle panel 
bears out the expected particle mass scaling when V2 is plotted vs. KEj. For KEj > 1 GeV, clear 
splitting into a meson branch (lower V2) and a baryon branch (higher vq) occurs. However, both 
of these branches show rather good scaling separately. The right panel shows the resulting values 
obtained after scaling both V2 and KEj (ie the data in the middle panel) by n q . This latter scaling is 
an indication of the inherent quark-like degrees of freedom in the flowing matter. That is, we assert 
that the bulk of the elliptic flow develops in the pre-hadronization phase. 



3.0.1 Validation of the pre-hadronic development of elliptic flow 

Despite robust quark number scaling, one can still ask whether or not the bulk of the elliptic 
flow develops pre- or post-hadronization. This question can be readily addressed via careful study 
of the V2 for the <p meson. The <p meson is comprised of a strange (s) and an anti-strange (s) quark 

2 Validation of the hydrodynamic scaling prediction that V4 <>= v| is presented in Ref. It is a non-trivial prediction; 
the measured differential \>2(pt) show a linear pj dependence for low pj pions but va,(pj) is quadratic in pj. 
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1 1.05 1.1 1.15 IS 1 1.05 1.1 1.15 1.2 

m inv (GeV/c*) m inv (GeV/c*) 

Figure 5: (Color on line) An example of invariant mass distributions for K + K~ pairs identified in the 
PHENIX detector for = 1.6-2. 7 GeV/c and reaction centrality 20-60%: The left panel shows the distri- 
bution for signal plus combinatoric background. The right panel shows the distribution after background 
subtraction. 



and its interaction with nuclear matter is suppressed according to the Okubo-Zweig-Izuka (OZI) 
rules [43]. Therefore, the meson is expected to have a rather small hadronic cross section with 
non-strange hadrons (~ 9 mb) [ 44 , 45|, [f6|] and this leads to the relatively large mean free path X§ , 
when compared to the transverse size of the emitting system. 

Thus, if elliptic flow was established in a phase involving hadrons interacting with their stan- 
dard hadronic cross sections (post-hadronization), one would expect V2 for the to be significantly 
smaller than that for other hadrons (eg. p and %). An added bonus is that the relatively long lifetime 
(~ 45 fm/c) for the (p makes its decay inside the pre-hadronization medium unlikely and this makes 
for a cleaner measurement. For these reasons, flow measurements for the play a crucial role for 
the distinction between pre- and post-hadronic development of elliptic flow [47]. 

The elliptic flow measurements for the <p meson are made difficult due to combinatoric back- 
ground; for the case presented here, such a background results from K + K pairs which do not 
come from the decay of the (p. Figures |5ja and b show a typical example of invariant mass dis- 
tributions for K + K pairs obtained in Au+Au collisions at i/sJviv = 200 GeV. These preliminary 
data show results before and after combinatorial background subtraction. Here, it is important to 
stress that the combinatorial background also exhibits an azimuthal anisotropy (V2) which needs to 
be separated from that for the meson. This can be achieved via a technique recently developed 
by Borghini et al. Jig). The technique which is based on the study of V2 as a function of inavariant 
mass, allows for robust V2 extraction even for a rather small signal-to-background ratio. Needless 
to say, the distributions shown in Fig |5] allow a very reliable V2 measurement for the (j) meson. 

The left and right panels of Fig. || compares the unsealed and scaled results (respectively) for 
V2 vs. KEj for 7i,K,p and <p in 20-60% central Au+Au collisions at ^/snn = 200 GeV. These are 
recent preliminary results obtained by the PHENIX collaboration. The left panel clearly shows 
that, despite its mass which is comparable to that for the proton, V2(KEj) for the <p follows the flow 
pattern of the other light mesons (n and K) whose cross sections are not OZI suppressed. The right 
panel indicate the expected universal scaling after n q scaling. We therefore conclude that, when 
elliptic flow develops, the constituents of the flowing medium can not be ground state hadrons 
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Figure 6: (Color on line) Comparison of V2 vs. KEj for K,K,p and <j) mesons (left panel) detected in semi- 
central (20-60%) Au+Au collisions at ^/snn = 200 GeV. The right panel show the same data scaled by quark 
number. These preliminary data have been obtained by the PHENIX collaboration. 



interacting with their standard hadronic cross sections. Instead, they reflect a pre-hardonization 
state of the created high energy density matter that contain the prerequisite quantum numbers of 
the hadrons to be formed ie. the QGP. It is noteworthy that a similar investigation can be made via 
V2 measurements for "heavy" multi-strange baryons such as the Q. + Q.; they too, are expected to 
have a relatively small hadronic scattering cross section. Indeed, we have observed that the recent 
high-statistics preliminary STAR results for CI + CI V2 in Au+Au collisions at y/sfw = 200 GeV 
(presented at SQM 2006) pi], 50] do follow the universal V2/n q vs KEj/n q scaling (see Fig Eh. 



3.1 Validation of charm quark diffusion and flow 

As indicated earlier, perfect fluid hydrodynamics stipulates that all constituents of the QGP 
should experience the same velocity field. This requirement places severe constraints on the flow 
of heavy quarks and their attendant hadrons because their relaxation time Xr is much larger than 



that for light quarks [49]; 



M 

Xr ~ ^r,, 

where M is the mass of the heavy quark, T is the temperature and Y s = T] / sT is an estimate of 
the light quark relaxation time; r\ and s are the viscosity and entropy density respectively. For the 
charm quark (M ~ 1.4 GeV) at a temperature of 165 MeV, this translates into an ~ 8 fold increase 
compared to T s . Heavy quarks are also produced with a power-law or non-thermal transverse 
momentum spectrum. Thus, one might expect the elliptic flow for heavy quarks to be much smaller 
than that for light quarks (for reasonable system lifetimes), unless in-medium interactions are strong 
enough to enforce many interactions with sufficient momentum transfer to thermalize them. 
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Figure 7: (Color on line) (a) Evolution of charm quark transverse momentum spectrum with different values 
of the diffusion coefficient D c as indicated. The thermal spectrum for T = 165 MeV and the initial transverse 
momentum spectrum given by leading order perturbation theory (LO pQCD) are indicated, (b) Evolution of 



charm quark v%(pt) with D c . These figures are taken from Ref. [49| 



The heavy quark diffusion coefficient D/ iq , is related to Tr p9^]; 

M 

Thus, Dj iq ultimately controls the extent to which the initial power-law spectrum approaches the 
thermal spectrum and the extent to which the heavy quark will follow the underlying flow of the 
medium. This is illustrated for the charm quark in Fig. \j\ [ p9| ] where the evolution of the pj 
distribution and the differential v%{pr), are shown for different values of the charm quark diffusion 
coefficient D c . The figure clearly indicates that measurements of charm quark differential elliptic 
flow V2(pr) can serve to constrain D c . 

One method of accessing charm v-i is to measure the the elliptic flow of "non-photonic" elec- 
trons that originate primarily from the semi-leptonic decays of D and B mesons. Such a measure- 



ment has been made recently and found to yield significant anisotropy for D mesons []14[]. 

Figures ||a and b show differential V2{pr) results obtained for non-photonic electrons and 
D mesons respectively. The D meson results reflect estimates for pj < 2.0 GeV/c, obtained via 
detailed simulations constrained by the differential measurement V2(pt) for non-photonic electrons 



Q14p . For these simulations, the pj spectrum was constrained by data and all non-photonic electrons 
were assumed to come from D meson decay [[14]]. Estimates were obtained for three separate shape 
assumptions for Vzipj) for the D mesons. The result from each is show in |8|b. 

These results indicate clear evidence that the D meson and consequently the charm quark, do 
flow with a large v<i- This suggest that the in-medium interactions are strong and perhaps frequent 
enough to thermalize the charm quark. If this is indeed the case, then V2(KEt) for the D meson 
should show the same universal scaling presented earlier for other hadrons. 
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Figure 8: (Color on line) (a) v>2 vs. for non-photonic electrons. The solid and dashed lines compare the 
results from a calculation which includes/excludes charm flow [|o|]. (b) V2 vs. pj for D mesons derived from 
the non-photonic electron measurements for pj < 2.0 GeV/c. Results are shown for several assumptions for 
the shape of the pj dependence of the D meson Vi. The results for pions are shown for comparison. 



Figures ||a and b compare the unsealed and scaled results (respectively) for V2 vs. KEj for 
7l,K,p and D mesons measured in minimum bias Au+Au collisions. The assumed shape for the 
transverse momentum dependence of the D meson V2 is that for the proton cf. Fig. |8]b. We note here 
that this assumption is compatible with a large array of differential flow measurements for different 
"heavy" particle species (see for example Refs. [11, [13|, 31, 30]). Both panels (a and b) of Fig. 
|9]show robust scaling; they confirm that the D mesons (for pj < 2.0 GeV/c) and their associated 
charm quarks flow with a fluid velocity field identical to that of the other constituents. This means 
that the medium responds as a thermalized fluid and the transport mean free path is small. 

The filled triangles in Fig. ^ show results from an approximate interpolation of the model 
estimates for the charm quark p9| ] presented in Fig. [7]b. In Fig. ^ these estimates are shown 
only for pj < 2.0 GeV/c ie. the available range of the D meson data 3 ; they are close to the ones 
obtained from the LO-QCD charm quark spectrum for a calculation in which D c ~ 3/2nT (see. 
Fig. 0b). We note here that the calculations are for the charm quark and do not include the effects 
of hadronization/fragmentation leading to D meson production. Despite the many caveats of the 
model [^], Fig. || shows that V2 for the charm quark can be made to follow universal scaling for an 
appropriate choice of the coefficient D c . It is noteworthy that an even smaller value D c ~ 1.5 /2kT 
is found if the drag on the charm quark is made to depend on the velocity ie. dp/dt °< v. 



4. Transport properties of the high energy density fluid 

The observed scaling of elliptic flow provide important constraints for the transport properties 



3 The indicated scaling is observed to break for charm quarks with pj > 2.0 GeV/c. Thus, high pj measurements 
are critical to the development of further constraints. 
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Figure 9: (Color on line) (a) v% vs. KEj. (b) V2/n q vs. KEj jn q . Results are shown for several particle 
species produced in minimum bias Au+Au collisions at ^snn = 200 GeV. Interpolated results from a model 
calculation |W9p for charm quarks are also indicated. 



of the high energy density QGP fluid. In what follows, we apply several of these constraints to 
make estimates for a number transport coefficients. 



4.1 Time scale for the development of elliptic flow 

As discussed earlier, elliptic flow develops gradually in the evolving system on a time scale % 
characterized by its transverse size R and the sound speed c s : 



R 

T ~ — . 

C s 



(4.1) 



Values for R have been reported as a function of centrality in Ref. Q25p for Au+Au collisions at 
y/SNN = 60 — 200 GeV. They indicate little change over the centrality range of interest here, so 
we use an averaged value (R) ~ 1.9. The value (c s ) = 0.35 ±0.05 determined via eccentricity 
scaling was presented earlier. These values give (t) ~ (R) / (c s ) ~ 5.4 fm, ie. the mean time re- 
quired to "translate" the initial spatial anisotropy of the high energy density fluid into a momentum 
anisotropy. 

4.2 The viscosity to entropy ratio r\/s 

Validation of the predictions from perfect fluid hydrodynamics for the scaling of the elliptic 
flow demands a low shear viscosity to entropy ratio (T] /s) for the high energy density matter created. 
On the one hand, this is very well supported by the observation that experimental V2 values are in 
fairly good agreement with simulated values based upon ideal hydrodynamics [ flOj |4T| , 21, p2| ]. 
On the hand the predicted suppression of flow, due to shear viscosity, by weak-coupling transport 
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calculations is not observed. However, Tf/s cannot be arbitrarily small because quantum mechanics 
limits the size of cross sections via unitarity. An absolute lower bound is 1 /47i, reached in the strong 
coupling limit of certain gauge theories. This bound has been recently conjectured to hold for all 



substances [53] 



The observed universal scaling of V2 can be used to constrain a reliable estimate for rj/s 



[g4|,|55I]: 

- ~ TX f c s , (4.2) 
s 

where T is the temperature, Xf is the mean free path and c s is the sound speed in the matter. The 
temperature T = 165 ± 3 MeV was obtained via a hydrodynamically inspired fit to the scaled data 
shown in Fig. 0. The functional form used for the fit is I\ (w) /Iq(w), where w = KEj /2T and I\ (w) 



and Io(w) are Bessel functions [p4]. For c s , we use the estimate c s = 0.35 ±0.05, given earlier 
and in Refs. [28] and The mean free path estimate Xf = 0.3 ±0.03/wj, was obtained from an 



on-shell transport model simulation of the gluon evolution (in space and time) in Au+Au collisions 



at 200 GeV [p6p. This parton cascade model includes pQCD 2^2 and 2 <-> 3 scatterings. These 
values for Xf, T and c s lead to the estimate tj /s ~ 0.09 ± 0.015 with an estimated systematic error 
of +0. 1 primarily due to uncertainties in Xf. 

This estimated value for rj/s is in good agreement with the experimentally based estimates 



of Teaney and Gavin Q57|, [58]] and the theoretical estimates of Gyulassy and Shuryak \p2[ p9j ]. 
However, it is much lower than the value obtained from weak coupling QCD [ jo"!] ] or hadronic 
computations [|60|]. This low value for r]/s has been interpreted as evidence that the QGP is more 



strongly coupled than previously thought [40, 52] 



4.3 The bulk viscosity 



Two component models (see for example Ref. Q51|]) are often used to estimate the bulk viscos- 
ity g. Assuming such a model holds here, the bulk viscosity can be expressed in terms of the shear 
viscosity and the sound speed as: 

? = ^(l-3c?) 2 . (4.3) 

Using the value (c s ) = 0.35, we find the relation (g) ~ 0.7T]. Thus, the implication is that the 
time averaged bulk viscosity of the high energy density matter is also relatively small but non-zero. 
Here, it is important to stress that for higher temperatures (eg. T ~ 2T C ), lattice QCD calculations 
give c 2 K ~ 1/3 which leads to g ~ 0. This implies that much of the bulk viscosity develops at the 
lower temperatures close to T c . 

4.4 The sound attenuation length 

With T and r\/s in hand, the sound attenuation length Y s , can be evaluated as: 

4 r? 

r, = ^. (4.4) 
3 sT 

This gives Y s ~ 0.16 fm and f]/Ts ~ 0. 12 fm. These values suggests a rather short relaxation time 
for light quarks. 
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4.5 The Reynolds R e and Knudsen K n numbers 

The Reynolds R e and Knudsen K n numbers play an important role in hydrodynamic consider- 
ations; 

(4.5) 

n 

Here £ is the energy density. A large Reynolds number indicates that viscous forces are not impor- 
tant to the flow. That is, the smallest scales of fluid motion are undamped. 
The Knudsen number is given as: 

K n = |; (4.6) 
K 

its reciprocal gives a measure of the number of collisions per particle. Thus, Local thermal equi- 
librium is expected if K~ l » 1. 

Since the shear viscosity v\ can be expressed as: 

T] = eXc s , (4.7) 

and flow is transonic (ie. Mach number M n ~ 1), 

R e ~^r. (4.8) 

Using the previously obtained values for X and R we obtain K ~ 6. Such a large number favors 
rapid local thermalization. 

4.6 Diffusion coefficient D c for the charm quark 



The diffusion coefficient D c , for the charm quark is related to its relaxation time Xr [19]; 



M 

Tr = -D c . (4.9) 

As discussed earlier, it controls the extent to which the initial power-law spectrum for the charm 
quark approaches the thermal distribution and the extent to which it follows the underlying flow of 
the medium. 

The agreement between theory and experiment shown in Fig. |9|, was achieved for 

D °~if (410) 

This gives the values D c ~ 0.6 fm and Tr ~ 5 fm. Thus, D c is approximately a factor of five times 
larger than the hydrodynamic scale rj/sT but Xr is similar to the time for the development of elliptic 
flow. 

4.6.1 Does the value of rj/s signal the QCD critical end point? 

Our earlier estimate of the viscosity to entropy ratio is indicated by the filled hexagon is shown 



in Fig. 10. This figure shows rj/s vs (T — T c ) /T c for molecular, atomic and nuclear matter. The 



data for He, N2 and H2O were obtained for their respective critical pressure. They are taken from 



Ref. [66] and the references therein. The calculated results shown for the meson-gas (for T < T c ) 
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Figure 10: (Color online) rj /s vs (T — T c ) /T c for several substances as indicated. The calculated values 
for the meson-gas have an associated error of ~ 50% [p2[], The lattice QCD value T c = 170 MeV J63| | is 
assumed for nuclear matter. The lines are drawn to guide the eye. 



are obtained from chiral perturbation theory with free cross sections [|62|]. Those for the QGP (ie. 
for T > T c ) are from lattice QCD simulations [J67Q. The value T c ~ 170 MeV is taken from lattice 
QCD calculations [|3|]. 

Figure [K] illustrates the observation that for atomic and molecular substances, the ratio j]/s 
exhibits a minimum of comparable depth for isobars passing in the vicinity of the liquid-gas critical 



point [ |53[ p6| |. When an isobar passes through the critical point (as shown in Fig. |10|), the minimum 
forms a cusp at T c ; when it passes below the critical point, the minimum is found at a temperature 
below T c (liquid side) but is accompanied by a discontinuous change across the phase transition. 
For an isobar passing above the critical point, a less pronounced minimum is found at a value 
slightly above T c . The value r\/s is smallest in the vicinity of T c because this corresponds to the 



most difficult condition for the transport of momentum [ 66p. 

Given these observations, one expects a broad range of trajectories in the {T,hb) or (T,/j.b) 
plane for nuclear matter, to show v\ / s minima with a possible cusp at the critical point. The exact 
location of this point is of course not known, and only coarse estimates of where it might lie are 
available. The open triangles in the figure show calculated values for r\/s along the = 0, ns = 
trajectory. For T <T C the t]/s values for the meson-gas show an increase for decreasing values of 
T. For T greater than T c , the lattice results [^] indicate an increase of v\ /s with T, albeit with large 
error bars. These trends suggest a minimum for r]/s in the vicinity of T c . This minimum is rather 
close to the absolute lower bound of f]/s = \/4n. We therefore speculate that it is compatible 
with the minimum expected if the hot and dense QCD matter produced in RHIC collisions follow 
decay trajectories which are close to the QCD critical end point (CEP). Such trajectories could 
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be readily followed if the CEP acts as an attractor for thermodynamic trajectories of the decaying 
matter @,^5j]. 

5. Summary 

In summary, elliptic flow measurements at RHIC provide compelling evidence for the 
production and rapid thermalization (less than 1 fm/c after impact) of a new state of matter having 
an energy density well in excess of the critical value required for de-confinement i.e the QGR 
Validation of the predictions from perfect fluid hydrodynamics for the scaling of the elliptic flow, 
in conjunction with the universal scaling of baryons and mesons (via quark number), serve as 
constraints for the extraction of several transport coefficients for this matter. These coefficients 
give important insights on the transport properties and the decay dynamics of the QGP and 
suggest a plasma having essentially perfect liquid-like properties. 
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